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Abstract 

Phase transitions in open quantum systems, which are associated 
with the formation of collective states of a large width and of trapped 
states with rather small widths, are related to exceptional points of 
the Hamiltonian. Exceptional points are the singularities of the spec- 
trum and eigenfunctions, when they are considered as functions of a 
coupling parameter. In the present paper this parameter is the cou- 
pling strength to the continuum. It is shown that the positions of the 
exceptional points (their accumulation point in the thermodynamical 
limit) depend on the particular type and energy dependence of the 
coupling to the continuum in the same way as the transition point of 
the corresponding phase transition. 
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1 Introduction 



Recently, mechanisms of restructuring of quantum systems are discussed with 
renewed interest HH - • Conditions for the formation of collective states on 
the one hand and of quantum chaos on the other are being studied. Among 
the oldest examples of the former is the schematic model explaining in simple 
terms the origin of giant dipole resonances in nuclei |p. The Hamiltonian is 
of the form H = Hq + Hi = Hq + DD^ where Hq is the Hamiltonian of the 
unperturbed system and DD^ describes the factorized dipole-dipole residual 
interaction. The rank of Hq is N where N is the number of unperturbed 
states considered in a certain energy interval, while the rank of DD^ is 1. 
If the modulus of the average matrix elements of the vector D is sufficiently 
large (Z) ^ do, where D is the average matrix element of D and do is the 
mean level distance of the eigenvalues of Hq), one of the eigenstates of H is 
shifted considerably in energy. Since its eigenvector has contributions from 
almost all basis states, it is a collective state. This type of collectivity is 
called internal collectivity . 

Besides the internal collectivity, there exists the so-called external collec- 
tivity of resonance states It appears at high level density since the discrete 
states are in general embedded in the continuum and coupled to each other 
via the continuum. The Hamiltonian of the open quantum system is, to a 
good approximation, given by Hq — iVV^ where VV'^ contains the coupling 
matrix elements between the discrete states j and the decay channels c. 
The rank of VV'^ is equal to the number K of open decay channels. It is 
usually much smaller than the rank N oi Hq. If the matrix elements of VV'^ 
are sufficiently large {V 3> (io), K collective states are formed. They are 
distinguished from the N — K non-collective states by their large widths or 
short lifetimes [0] - 101 • The wavefunctions of these K collective resonance 
states are again characterized by a large number of components of the basis 
states. 

There is a basic difference between internal and external collectivity. In- 
ternal collectivity deals with a self-adjoint Hamiltonian, where the coupling 
parameter is real, while external collectivity is described by a purely imag- 
inary coupling strength. The particular type of collectivity is therefore ex- 
pressed by a real spectrum in the first case, whereas in the second case the 
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dissipative character of the open system is reflected by imaginary parts of 
the eigenenergies, the physical widths. Formally, the formation of a collective 
state is, however, the same in both cases. It is expressed by the mixing of 
the eigenf unctions of the full Hamiltonian with respect to the eigenf unctions 
oiHo. 

The case rank(ifo) = rank(i7i) is not necessarily connected with the ap- 
pearance of collective states. It has been widely discussed in the hterature 
[P]. Conditions have been given for the occurrence of chaotic behaviour and 
of some sort of phase transitions |T^. The mechanism of such restructuring 



and possible chaotic behaviour in the transitional region (A ~ Ac) have been 
explained in terms of the exceptional points (EPs) of the problem Hq + AHi. 

The EPs are the only singularities in the complex A— plane of the eigen- 
values Ek{A). Their positions are fixed by the choice of Hq and Hi only. As 
a consequence, the distribution of the EPs is characteristic for any particular 
Hamiltonian of the form Hq + AHi, once Hq and Hi are given. 

It is therefore natural to discuss the EPs as they determine important 
properties of the spectrum such as the statistical properties, the type and lo- 
cations of avoided level crossings, the softness of the spectrum and the ranges 
of A— values, where special features of the spectrum occur. In particular, the 
distribution and frequency of the EPs give indications about transitional re- 
gions and possible occurrence of chaotic behaviour. 

The studies of the EPs have, in the papers quoted, been dealing with 
closed systems, where the influence of decay channels in the continuum has 
not been considered, that is H is selfadjoint and the eigenstates are dis- 
crete. As stated above we are often faced with dissipative systems, where 
the coupling to the continuum and decay channels is expressed by complex 
eigenenergies which give rise to finite life times of resonance states. The 
complex eigenvalues £k = Ek — i/2Tk are interpreted as resonance states at 
the energies Ek with the decay widths which are inverse life times. 

While there is a great mathematical similarity in the treatment of the 
respective self-adjoint and dissipative Hamiltonians, the physical findings 
deserve particular attention. Of special interest will be the circumstances. 
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under which the formation of collective states can be understood as a phase 
transition. We generalize the results obtained in to open quantum me- 
chanical systems and expand the investigations of ||^ by taking into account 
the effects and properties of the EPs. In all cases, we will restrict ourselves to 
the one-channel case in which the phase transition as well as the collectivity 
of one of the states are well pronounced. 

Generally, a phase transition is a substantial restructuring of the system 
taking place at a finite (critical) value Ac of a certain control parameter A. 
This restructuring is a collective phenomenon, ranging over all scales inher- 
ent in the system. The nature of the reorganization process is characterized 
by the behaviour of an order parameter being a function of A. In our case of 
an open quantum system the order parameter is Tq/N where Tq is the width 
of the collective mode (in the one channel case) and is the total number 
of states. Its first derivative with respect to A shows a finite discontinuity 
at the critical value Ac corresponding to a second-order phase transition. All 
characteristic features of the phase transition emerge already at finite values 
of despite the fact that the strict thermodynamical definition is possible 
only for A^ — oo. Precise conditions for its occurence are derived in [0. 

The phase transition is connected with the appearance of a state whose 
external collectivity is of a global nature. Its wave function carries contribu- 
tions from a// eigenstates of Hq. This global collectivity must be distinguished 
from the local (external) collectivity which appears when the resonance states 
are coupled strongly to the continuum but the conditions for a phase tran- 
sition are not fulfilled. The wave function of this collective state carries 
contributions only from a restricted number of basis states, whose unper- 
turbed energies are overlapped by the width of the collective state ||^. 

In Section 2 the basics of the EPs and specific features relating to open 
quantum systems are reviewed. A simple example using two resonances cou- 
pled to one decay channel illustrates the connection between the main results 
known from the study of open quantum mechanical systems and the aspect 
of the EPs. In Section 3 the restructuring of the system is put into context 
with several systems of a different nature. As a beginning, the simplest case 
of a picket fence model is considered. Characteristic patterns of the EPs 
relating to phase transitions are established. Next, the influence of random 
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perturbations within the system is investigated. Finally we discuss the effect 
of a more general level density distribution of the unperturbed levels and that 
of the coupling vector V. A short summary and outlook is given in Section 4. 



2 Exceptional points and features of open 
quantum systems 

The N X N matrix problem of the form Hq + AHi has eigenvalues Ek{A), k = 
1, . . . ,N which are obtained by the secular equation 

det{Ho + AHi -E) = 0. (1) 

We here assume that Ho and Hi are real and symmetric. The EPs are 
characterized by the coalescence of any two pairs of eigenvalues. They are 
fixed by a particular choice of Hq and Hi. If we exclude a genuine degen- 
eracy of eigenvalues for real values A, such coalescence will happen only for 
complex values of the coupling parameter and hence for complex eigenvalues. 
Accordingly, the EPs are determined by the simultaneous solutions of Eq. (|^) 
and of 

^det{Ho + AHi-E) = (2) 

where the second equation ensures that two eigenvalues coincide. Eqs.(||) and 
(0) are polynomials in E and A of order N and A^ — 1, respectively. They can 
be combined into a single polynomial in A of order A^(A^ — 1), the resultant 
of Eq.(|T]), by eliminating the variable E. The resultant has real coefficients, 
hence the EPs occur in A^(A^ — l)/2 complex conjugate pairs. They are the 
only singularitites which the eigenvalues Ek{A) can have as functions of A. 
In fact, they are the square root branch points of one analytic function which 
has N Riemann sheets, where the values on each sheet are just the real eigen- 
values Ek{A) for real values A. All eigenvalues are therefore analytically 
connected with each other via the square root branch points in the complex 
A— plane. 

The behaviour of the eigenfunctions when continued to the EPs deserves 
special interest. Since they are lying in the complex A-plane, the operator 
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Hq + AHi is no longer self-adjoint at A = Aep with Aep denoting a (complex) 
EP. Therefore we cannot expect two linearly independent eigenvectors even 
though two (complex) eigenvalues Ei and coalesce. The eigenbasis of Hq + 
A^pHi is no longer orthogonal as it is the case for selfadjoint Hamiltonians. In 
fact, in contrast to the self-adjoint case, where a twofold degeneracy always 
implies a two-dimensional eigenspace, an EP is characterized by the fact, 
that the rank of the associated matrix Hq + Agp-ffi drops by one. In other 
words, there is no two-dimensional subspace associated with the coalescence 
of two eigenvalues. We rather encounter a confluence of the two eigenvectors 
ipi{A.) and ipk{-A) at A = Aep- Moreover, for a general complex value A not 
coinciding with an EP, we may choose a biorthogonal system such that 

(MmkiA)) = (3) 

where ipi{A.) and V'i(^) ^ire the left and the right eigenvectors, respectively. 
We note, that the transformation O diagonalizing a symmetric complex 
Hamiltonian H is complex orthogonal {00^ = 0^0 = 1) but not uni- 
tary. Therefore, the scalar product of the two left or right eigenvectors obeys 
(?/'fc(A)|?/'fc(A)) = (^/'fc(A)|-?/'fc(A)) > 1 for any k. A problem occurs when the 
normalization condition given by Eq.(^ is considered for A — > Aep: with 
Ei{Aep) = -Efc(AEp) it is ^/'z(Aep) = ipk{AEp)- This means, in view of Eq.(|), 
that ipk{AEp) cannot be normalized at A = Aep, since now the orthogonal- 
ity conflicts with the normalization requirement. Usually, Eq.(|^) is globally 
enforced; as a consequence, the two states i^i^A) and ipkiA) not only coin- 
cide for A Aep but they blow up, that is limA^AEp(V^fe(^)lV^fc(^)) ~^ [|]. 

The physical significance of the EPs lies in their relation to avoided level 
crossing. In particular, in a number of examples it has been demonstrated 
that in the region of the real values of A, where a high density of EPs occurs, 
the typical statistical characteristics of the spectrum ascribed to quantum 



chaos prevail P, ITT 



Turning to open quantum systems 

H = Ho- iAVV^ (4) 

we are faced with a Hamiltonian which has in general complex eigenvalues 
Sk = Ek — i/'iTk as A is complex, 

A = \e'^ . (5) 
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(The factor i in front of the couphng term A in Eq. (^) is used tradition- 
ally in this context). The question then arises, whether and under which 
conditions two (or more) eigenvalues coincide in the complex energy plane 
and how such a crossing depends on the coupling strength A. The answer to 
these questions depends on the manner by which the EPs are fixed by the 
operators Hq and V. 

For illustration let us consider the simple example of two resonances, 
which are coupled to one open decay channel. The Hamiltonian matrix for 
this system can be written in the eigenbasis of Hq as 

/ ei \ . . / cos^ uj cos uj sin u; \ 

H = { n - • -2 

V €2 / V COS uj sm io sm lo ) 

= H^-ikVV\ (6) 

where we use = VV'^ . The relative coupling strength of the two resonance 
states to the continuum is determined by the angle u of the vector V = 
(coscj, sincj). The simplicity of the model provides an analytic expression 
for the two EPs, viz. 

Aep = ^(€2 - ei)e±2- (7) 

which are the zeros of the square root in the expression for the eigenenergies, 
which read 

gi.2 = '^ + '^^-'^ ± y^,^ _ ,^)2 _ 2^A(e, - 62) cos 2a. + (zA)^. (8) 

The example nicely demonstrates that the EPs depend only on Hq and Hi. 
Under variation of the angle u the EPs are moving on a circle with radius 
62 — ei in the complex A— plane. What determines the EPs are the energies 
of the unperturbed states (that is Hq) and their relative coupling (that is 
Hi), and not the phase of the coupling A = Ae**^. 

Yet, the complex eigenvalues £^1,2 obviously do depend on A, that is on 
its modulus A and on the angle ip. Variation of A invokes trajectories in the 
complex energy plane. The contours of the trajectories depend therefore on 
both, the phase ip and the relative coupling controlled by u. For u = 45° and 
f = V^EP = 0° the two trajectories will cross one EP at A = Aep = £2 — ei. At 
this critical value of the coupling strength the two eigenvalues coalesce. The 
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dimension of the eigenspace of the Hamiltonian is reduced to one including 
the consequences for the eigenvectors mentioned above. We stress, however, 
that the crossing of trajectories should not be confused with the familiar 
degeneracy for self-adjoint Hamiltonians. 

For 99 7^ (p-Ep the two trajectories repel each other within a certain finite 
distance in the complex plane @]. This is a generalization of avoided level 
crossing for real eigenvalues of a hermitian Hamilton operator. Only if if is 
properly tuned to the value of u can the two energy trajectories genuinely 
cross each other. A nice illustration of this simple example for different val- 
ues of u and ip can be found in . 

What becomes obvious in the 2x2 matrix model can be generalized to an 
arbitrary N x N situation, i.e. resonance states coupled to one common 
decay channel. With increasing A one eigenvalue trajectory always drifts fur- 
ther into the complex plane, while the others are bending back towards the 
real axis after they have repelled with the collective state whose width always 
increases. This happens irrespective of ip. Physically it means that one of 
the resonance states takes almost all of the transition strength by trapping 
the others which then become long-lived 0, |^ . 

This is understood by the rank one of the coupling matrix. At large values 
of A the first part Hq of the Hamiltonian is a small perturbation. Therefore, 
the total matrix is essentially turned into an operator of rank 1. In other 
words, there is only one nonzero eigenvalue, and the widths of the zero eigen- 
values have to vanish. In general, with A^ resonance states and K < N open 
decay channels, there appear K fast decaying states and N — K states which 
are virtually stable. Various examples of the many channel - many resonance 
case can be found in - 0. 

The low rank of Hi has a drastic effect on the total number of EPs. If K 
is the rank of Hi , the number of EPs is K{2N — K — 1). The important point 
is the linear behaviour in A^; only when K attains the order of magnitude 
of A^ is the quadratic behaviour retrieved. This finding is significant in that 
it indicates that chaotic behaviour cannot be generated by a low rank of Hi, 
there is simply an insufficient number of avoided level crossings. This is in 
line with an analysis of the level statistics of a Poisson ensemble coupled to 
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a continuum by a Gaussian coupling vector 0, ^, where, for large coupling 
strength, chaotic features of the system are the more pronounced the larger 
the number K of decay channels. 



3 Results 

We turn our attention to the critical region, where a restructuring of the sys- 
tem occurs. We address the relationship between the potential occurrence of 
a phase transition which is associated with the formation of a globally col- 
lectice state and the distribution of the EPs. Since the essential properties of 
the Hamiltonian must be reflected in the distribution of the EPs, we expect 
this relationship to exist. 



3.1 Simple example: picket fence 

We consider the simple model with 
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and the N x N matrix of rank one 



(I 



Vi 



(9) 



(10) 



The Hamiltonian describes a picket fence spectrum, where all states are cou- 
pled equally to one decay channel. Such a system shows a phase transition at 
A = I/tt 0. Results for finite N can be obtained easily by numerical means. 
In the limit N ^ oo the zeros of det{Ho—iAVV'^—S) and sin{'n'S)—iA cos(7r£) 
coincide [^. One finds an accumulation point of the EPs, which emerges in 
the limit and is found from the zeros of the resultant, which is obtained by 
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eliminating the variable 8 from the set 



sin(7r£) +i7rAcos(7r£) = (11) 
cos(7r^) -i7rAsin(7r^) = 0. (12) 

These two equations are equivalent to Eqs. (|I]), (0). Their simultaneous 
solutions, the zeros of the resultant (denoted by Rsl), are given by 



Rsl(A) = VI + (^7rA)2. (13) 

Obviously, this is no longer a polynomial. The important point of our finding 
is that all roots which occur in the complex A— plane for finite converge 
to A = A = il/vr in the limit oo. Negative values of A would lead 

to negative decay widths which have no physical meaning. In the following, 
we restrict ourselves to positive A and 0° < < QO'', since all relations are 
symmetric with respect to the replacement A — > —A and A — > A* . Note also 
that the limit point, being a point of accumulation, is no longer a square root 
branch point for the energy spectrum but rather a logarithmic branch point. 
In fact, the analytic behaviour of the (infinitely many) energy levels ^fc(A) 
in the vicinity of A = I/tt is found by solving the secular equation explicitly 
for £. The expression reads 

f(A) = -atanHAx) = -l„— (14) 

which clearly reveals the logarithmic branch points at A = l/vr. Moreover, 
for Avr < 1 we read off 

Sk = k-ik + 0{2) (15) 



while for Avr > 1 



Sk = k -\ 1 — -In 

2 2tx 



1 - Att 



1 + Avr 



with k integer. In both cases, the remaining terms, denoted by O, are purely 
imaginary if = 0°. These results were obtained in [^. 

For finite N the resultant relating to the Hamiltonian Hq — iAVV'^ be- 
comes a polynomial of order — 1 in A^. The complex roots which are the 
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EPs are therefore not only symmetric with respect to the real axis but a so- 
lution Aep implies also the solution — Aep- As a consequence, a solution has 
to occur on the real A— axis for N even. This is a non-generic feature of the 
present model. For odd values of and (f = 0°, no eigenvalue is crossing an 
EP under variation of A. We use such cases in our numerical demonstrations. 

To illustrate how the S undergo level repulsions and how the repulsions 
are related to the crossing points £{Aep), we display in Fig.l the crossing 
points and the eigenvalue trajectories for different values of A^ and (f = 0°. 
Since the spectrum is symmetric with respect to positive and negative en- 
ergies, only the positive part near to the centre is drawn. For A = all 
trajectories begin at the unperturbed energies e^, move into the complex 
plane and then turn back (e^ 7^ 0) towards the real axis again. For large A^ 
the corresponding values are given by Eq. (|T^ (which is valid for A^ 00). 
With increasing A^ the imaginary parts of the turning and the crossing points 
increase, while they move nearer to each other. The trajectory of the collec- 
tive state (efc = 0) is moving on the imaginary axis towards larger imaginary 
values implying an increasing width P. We interpret the turning points as 
level repulsions of the collective state with the other levels. The descending 
slope of the envelope of the turning points and the crossing points is due to 
the finite size of the spectrum. 

In Fig. 2 we display the EPs in the A— plane for a few values of A^. The 
sets of the EPs belonging to the same A^ are connected by a solid line. The 
zooming in of the EPs towards the accumulation point Ac = I/tt for increas- 
ing A^ is clearly discernible and enhanced in the insert of Fig.2. The larger A^ 
the larger is the density of the EPs near to the real axis. Additional points 
come in further away with each additional step of A^. They are typical edge 
effects and correspond to the turning points of the energy trajectories at 
the outer edges of the spectrum. These points which emerge further away 
from the real axis quickly move in to get near to the others while new further 
points come up for the following steps of A^. Eventually, in the limit N ^ 00, 
all points coalesce at the accumulation point. 

In this way, variation of the coupling strength A invokes a certain trajec- 
tory of the system in the complex A— plane. If it hits the accumulation point 
Ac (or the high-density regions of EPs at finite A^), all complex eigenvalues 
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coalesce (nearly coalesce) at their crossing point. This leads to a sharp tran- 
sition from a regime with resonance states to a regime with one collective 
state and — 1 trapped states. 

We note that at the accumulation point A = A^, which can occur only in 
the limit A^ —>■ oo, all infinitely many eigenstates collaps into one. It is the 
global collective state. It retains its characteristics also for A > Ac, when the 
other states (the trapped states) re-emerge. This particular case underlines 
the crucial connection between the critical point of a phase transition and 
the EPs. 

In the picket-fence model the accumulation point Ac lies on the real axis. 
If for physical reasons, a value f ^ has to be chosen with Hq and Hi left 
unchanged, a variation of A will effect a trajectory in the complex A— plane 
that passes the high-density region (or accumulation point for infinite A^) at 
a certain distance. This affects the sharpness of the transition between the 
two regimes. Note that this type of softening of the phase transition is dif- 
ferent in nature from the one caused by a finite value of A^. It persists in the 
limit N oo. This underlines that (f as well as the distribution of the EPs 
determine, whether there is a phase transition in the strict thermodynamical 
sense, which is associated with the sudden formation of a globally collective 
state in the system. 

As an aside we realize that the schematic model {(f = 90°) cannot have 
any signatures of a sharp phase transition, when the real value of A is varied 
from zero to large values. The accumulation point is then relatively far away 
on the imaginary axis of the coupling strength. Results for a variety of angles 
ip but the same Ho and Hi have been investigated in 0. 

Previously, a criterium for a phase transition, viz. 

1 ^ 

B = ^TW-T E (V^'cl^^fc) (17) 

-^^^ k=-N 

has been introduced in 0. Excluding the non-generic case of an acciden- 
tal crossing of two (or more) eigenvalues at an EP, B is an indicator of the 
sharpness of the transition. On the one hand, it is the larger the smaller the 
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distance between the complex eigenvalues and their crossing points. This 
distance is determined by the angle ip. On the other hand, B ^ 1 only if 
all complex eigenvalues reach their minimum distance simultaneously at the 
same value A = i.e. if the EPs accumulate. In this case, B shows a 
pronounced maximum as a function of A around A = Ag. 

In Fig. 3 a few cases including the picket fence model are illustrated. It 
also demonstrates a situation where a local collective state is formed, which 
means that the phase transition is washed out completely. We return to this 
latter aspect at the end of Section 3.3. 

To summarize the findings of this subsection: the high density of EPs 
provides the mathematical mechanism for the restructuring of the system 
under variation of the modulus A of the coupling parameter A towards larger 
values. The sharpness of the transition which is evoked by the variation is 
determined by the distance at which the high density region (accumulation 
point for infinite N) of EPs is passed by the corresponding energy trajecto- 
ries. This distance is determined by the angle </?. 



3.2 Random change of the unperturbed spectrum 

The essential aspects of our findings remain unchanged if the unperturbed 

energies in Hq and/or the elements of the coupling matrix V deviate from 
the symmetrical form used in the previous Section. We address the question, 
whether and to what extent such disturbance can change the basic pattern, 
that is the formation of a region of high density of EPs or even of an accu- 
mulation point in the large N limit. 

For this purpose, we define random perturbations in Hq by {Ho)k,k = 
— {N — 1 — 2k)/ 2 + Tfc where the are random numbers from a uniform dis- 
tribution in the interval [—0.1, 0.1]. In Fig.4 we have drawn the EPs for the 
perturbed and unperturbed picket fence in the complex A— plane for N — 19. 
As in Fig. 2, the unperturbed EPs are connected by a solid line. The per- 
turbed EPs do no longer lie on the smooth curve but they are scattered 
around it. 
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The quantity B of Eq. ( ]T7| ) indicates the accumulative behaviour of the 
EPs. It is illustrated in Fig.3 for the disturbed picket fence model as a func- 
tion of A for two different values of using = 0°. The excess beyond unity 
of B around A = l/vr increases with increasing N. We interpret this result 
as a strong indication for the disturbances to be washed out, the more so 
the larger N. In other words, for an increasing number of states, the average 
distribution of the EPs is zooming in to = Ac = l/vr. Hence, the EPs of 
the irregular system accumulate on the average. 

Furthermore, the precice form of the distribution from which the random 
changes in Hq and/or Hi are drawn is immaterial. In particular, choosing the 
eigenvalues of Hq from a Wigner (GOE) or a Poissonian distribution, does 
not alter our conclusions. The existence and the position of the accumulation 
point remains unaffected. 



3.3 Level density dependence 

In real physical systems, the level density, the number and the coupling 
strength of the decay channels are in general energy dependent. The former 
is usually a monotonically increasing function, whereas for instance in nuclear 
physics the continuum coupling strength decreases. This energy dependence 
will influence the distribution of the EPs. 

In it has been shown that it needs a proper tuning between the density 
dependence on the one hand and the coupling dependence on the other in or- 
der to guarantee a phase transition at a finite value of the coupling strength. 
Given the energy dependence of the level density p{E), a phase transition 
occurs at finite coupling strength, if and only if the energy dependence of 
the coupling vector is given by the inverse function p{E)~^. If the system 
obeys this condition on the average, a phase transition still occurs. The crit- 
ical point may be shifted with respect to the value of the ideal picket fence 
model. We talk about 'overcompensation' when the energy dependence of 
the coupling vector decreases (increases) at a lesser (faster) rate than that 
of the inverse behaviour of the level density. In this case, numerical results 
[0] led to the conjecture that the critical point is shifted to zero, i.e. Ac — 
if iV ^ oo. In the opposite case, which we denote as 'undercompensation'. 
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a global collective mode and a global reorganisation of the spectrum as a 
whole is absent. Now one obtains Xc ^ oo if N ^ oo. The occurence of the 
broad mode remains then a local phenomenon for all finite values of A [0. A 
collapse of the Hilbert space as in the picket-fence model does not occur. 



The following substantiates the conjecture about the position of Ac in the 
under- and overcompensated cases. To facilitate the discussion we restrict 
ourselves to a power behaviour of the energy dependence of the coupling and 
level density. The energy dependence of the coupling matrix elements {vkl"^ 
is assumed to be of the form \k\^. The unperturbed energies (eigenvalues of 
Hq) are defined as = sign(A;)|A;|*/^. The secular equation Eq. (|I]) can be 
written as 

k=-N/2 ^-^k A 
k=-N/2 



Using (|T8]), we obtain 



which corresponds to Eq. 



The EPs are simultaneous solutions of Eqs.(|l8l) and (|T9|). It is obvious 
from Eq.(p!9D that the energy values at the EP cannot be real. In fact, it 
is known [0, and also discussed in connection with Fig.l, that, for t = 2 
and r = (picket fence), the imaginary parts of the energies at the EP 
tend to infinity for — oo. Numerical evidence as well as the following 
consideration supports the conjecture that this holds for arbitrary values of 
r and t as long as t — r > 1. For this purpose we rewrite Eq.(|18]) for even A^ 



In the large A^ limit we replace the sum in Eq.(^0|) by an integral and £ by 
iz and obtain 



oo 



^'dk —2iz ^ ,^ ^ l + r„ 1+r 2\ 



■^'h TT^' = F373l^'-(U-— 2-— ;-.^) (21) 



I 

A' 
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The derivative with respect to z of the right hand side of Eg. (|2T|) , which is 
a hnear combination of two Hypergeometric Functions, corresponds to the 
large N hmit of Eq.(^) [E iz). It is straighforward to show, for instance 
by graphical means, that this derivative has no zero for finite values of z. In 
fact, the derivative vanishes only for z ^ oo irrespective of r and t {t > r + 1). 
We exploit this fact in Eq.(|18D in that we consider the limit of the left hand 
side for large imaginary values of S = E — ir/2. It is found to be proportional 
to r^^^~*"''^''*/r. This yields the generalised result that the overcompensated 
case occurs for 2(1 + r) > t, and the compensated and undercompensated 
cases for 2(1 + r) = t and 2(1 + r) < t, respectively 0]. In fact, the respec- 
tive limit of the left hand side of Eq.(^) is infinity, a finite constant and 
zero for F — > cx). Consequently, the corresponding values of A must be zero, 
finite and infinity, respectively. Since these values constitute simultaneous 
solutions of Eqs.(|TSp and (|T^), they are the values where the EP accumulate. 
Furthermore it follows, that the accumulation point is on the real axis for 
the compensated case {t = 2(1 + r)). 

To demonstrate further this result, the values of i? as a function of 
A = X {(f = 0) are drawn for the compensated case r = l,t = 4 in Fig.3. 
One clearly sees that the maximum of B around A = 2/7r increases with A^. 
As explained above, this means that, on the average, the minimum distance 
between the crossing point and its correponding eigenvalues decreases with 
increasing N. In the limit N oo, all eigenvalues will hit the crossing point. 
Thus, the compensated system behaves similar to the picket fence. Also in 
this case, the space of the eigenstates collapses at A = A^ = 2/7? into di- 
mension one, and a globally collective state is created. For comparison, the 
undercompensated case with r = and t = 4 is drawn in Fig.3 as an example 
for a system which does not undergo a phase transition. As explained above, 
Ac ^ cxD in this case and B has no maximum at a finite value of A. 



3.4 Phase behaviour of the wave functions 

In Section 2 some properties of the eigenfunctions have been pointed out, 
when a non self-adjoint Hamiltonian is considered. We resume the discus- 
sion in more detail here with particular emphasis on the phase of the wave 
functions when an EP is approached by the variation of the coupling pa- 
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rameter A. Surrounding one crossing point in the complex energy plane 
corresponds to a double loop around the EP in the A— plane, Eq. (J^), since 
the EPs are square root branch points, Eq. (^. We address in the following 
the effect of looping once around an EP in the A-plane. For demonstration 
we use the simple two dimensional model introduced in Section 2, Eq. (H). 

The two eigenf unctions, normalized according to Eq. @, can be parametrized 
by the complex angle 9: 

where the angle 9 is given by 

tan^ 9 = l^-l^-i^i-^^)+^^-o^^-_ 
Ci — £2 + (ei — 62) — iA cos 2iu 

The notation has been introduced in Eqs. @ and (H). From this expression 
we read off: 

(i) at an EP {Si = S2) we obtain tan^ 9 = —1 which implies |cos^| = 
I sin6'| = 00, that is the components of the wave functions blow up; 

(ii) when an EP is surrounded in the A— plane (which amounts to S1 — S2 

^2 — ^i); the tan^ 9 is changed into 1/ tan^ 9, which corresponds to the 
change 9 ^ 9 + n /2. 

This implies ipi ip2 and 1IJ2 — V'l- 

While it may not be obvious to implement such a contour in the complex 
A— plane in an actual physical experiment, there could be a possibility to 
achieve the same effect by a variation of the modulus of A for different values 
of the relative coupling given by u. In fact, using the setting Hq — i\VV\ 
then the EP lies for u > 45° just below and for u < 45° just above the point 
X = €2 — ci. If we compare the two situations while varying A, the two wave 
functions for A > 62 — ei differ in the same way as if the EP had been sur- 
rounded. 
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This can be made explicit by choosing an expression for tan 6' which is 
more convenient for this purpose, viz. 



— 2Asin2ix' 

tan^ = ^ . (24) 

ci — C2 + ei — ^2 — "^A cos 2u! 

The difference between the two values for u manifests itself in the differ- 
ence of the sign of the imaginary part of Si — S2, since different Riemann 
sheets have been approached. As a consequence, for large values of A the 
right hand side of Eq. (12^) tends towards tana; for the one case and towards 
— cotcij = tan(7r/2 + u;) for the other. 

We suggest that using electro-magnetic resonators may allow control of 
both, the global coupling A and the relative coupling between two resonances 
which is given hj u. If it should be possible to bring to interference the wave 
functions for the two different situations, that is for u > 45^ and for u < 45°, 
the different phases are expected to be observable. 

The situation described is reminiscent of Berry's phase |jl2|. However, we 



stress that in our discussion we are dealing with non self-adjoint operators, 
and an EP is therefore not to be confused with Berry's diabolic points. This 
difference was pointed out in In fact, a generic EP (i.e. when just 

two levels coalesce) gives rise to a double pole in the corresponding Green's 
function or S-matrix. The mechanism for this to occur is again related to 
the vanishing of the norm of the eigenvectors when A Aep. The double 
pole is an additional signature for an EP. It is different from the simple pole 
emerging in the Green's function (S-matrix), if a usual single resonance or an 
incidental degeneracy of two resonances occurs, say, owing to some symmetry. 
The latter situation is, however, an unlikely event as it needs the tuning of 
four parameters (six parameters in the case of no time reversal) to achieve 
it. In contrast, the existence of the EPs is an intrinsic mathematical feature. 



4 Summary and Outlook 

In this paper, we studied the relation between phase transitions, collective 
states and the distribution of EPs of open quantum systems. Since both the 
collectivity and the phase transition are well pronounced in the one-channel 
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case, wc restricted ourselves to an A^— dimensional quantum system coupled 
to one open decay channel. Our results are as follows: 

(i) A necessary condition for a phase transition to occur is that the EPs 
accumulate in the complex A— plane. This result has been shown analyti- 
cally for the picket-fence model and some generalisations, and is confirmed 
by numerical results using the quantity B. The quantity 5 is a measure 
for the average minimum distance of the eigenvalues to its crossing point 
in the complex energy plane. In the cases considered, B has a pronounced 
maximum at the critical value Ac at which the phase transition occurs. The 
increase of B with the number of states is a further indication that the 
EPs accumulate in the complex A— plane. 

(ii) The choice of the angle (f is dictated by the physical situation. Only 
if it is equal to the phase of A^, the system hits the accumulation point (or 
goes through the region of high density of EPs), when A is varied. Otherwise, 
the system may pass the vicinity of the accumulation point (or the region of 
high density of EPs) and a phase transition in the strict thermodynamical 
sense does not occur. In all cases considered, we find that (p has to be 0° 
for hitting the accumulation point Ac and for a genuine phase transition to 
emerge. 

(iii) If the system hits the accumulation point, the infinite-dimensional 
space collapses to a one-dimensional function space, since the infinitely many 
eigenfunctions are identical at A = Ac. This is in contrast to the function 
space of a usual A^— fold degeneracy, where N independent eigenfunctions 
occur. The eigenfunction of the collective state contains contributions of all 
basis functions of the unperturbed system Hq. It is a globally collective state 
which is created by the system as a whole. 

(iv) If the system docs not hit the accumulation point, the formation of 
the collective state is not related to a genuine phase transition. This is in 
particular the case for a system with a purely hermitian coupling matrix Hi 
{(p — 90°) which creates a state with large internal collectivity. 

(v) If the accumulation point is at Ac — oo, the formation of the collective 
state occurs successively and locally. Its wavefunction contains contributions 
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only from those resonance states which are overlapped by it. 

Exploiting the fact that the crossing points in the complex energy plane 
are square root branch points, we suggest a possibility for an experimental 
study of the level repulsion of the complex eigenvalues. For this, the relative 
coupling strength of the states to the continuum as well as the overall value 
of their coupling strength should be controllable. A signature of local reso- 
nance crossing is given by a particular change of the phase of the complex 
eigenfunctions. Work towards this aim is in progress. 
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Figure Captions 



Figure 1 

The trajectories of the eigenvalues for N — 15 and = 43 for increasing 
A e [0.001, 2] in steps of 0.001. The crossing points are marked by diamonds 
(A^ = 15) and crosses (A^ = 43). 

Figure 2 

The exceptional points in the complex A— plane for A^ = 15 (rhombs). A" = 19 
(plus signs), N — 27 (crosses) and A" = 43 (triangles). The inset is a magnifi- 
cation around the accumulation point Ac (black square) . The arrows indicate 
the changes of the EPs with increasing A^. 

Figure 3 

S as a function of A for different systems with (p — 0: ideal picket fence with 
A^ = 101 (dotted line), randomly perturbed picket fence with A" = 101, 1001 
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(solid lines), compensated case with r = 1, t = 4 and = 101, 1001 (dashed 
lines), undercompensated case with r = 0, t = 4 and iV = 101 (thick line). 
The two values Xc— I/tt and 2/7r, referring to A?" — > oo, are indicated by a 
vertical sohd line. 

Figure 4 

The EPs for A'^ = 19 for the ideal picket fence (rhombs) and the randomly 
perturbed picket fence (plus signs). Ac is denoted by a black square. 
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